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Abstract. We establish a unified model to explain Quasi-Periodic-Oscillation (QPO) observed from 
black hole and neutron star systems globally. This is based on the accreting systems thought to be 
damped harmonic oscillators with higher order nonlinearity. The model explains multiple properties 
parallelly independent of the nature of the compact object. It describes QPOs successfully for several 
compact sources. Based on it, we predict the spin frequency of the neutron star Sco X-1 and the 
specific angular momentum of black holes GRO J1655-40, GRS 1915H-105. 
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INTRODUCTION 

The origin of Quasi-Periodic-Oscillation (QPO) and its properties are still ill- 
understood. The observed QPO frequencies of compact objects are expected to be 
related to the spin parameter of the compact object itself. In certain neutron star systems 
a pair of QPO forms and QPO frequencies appear to be separated either by the order of 
the spin frequency of the neutron star, apparently for slow rotators e.g. 4U 1702-429, 
or by half of the spin frequency, for fast rotators e.g. 4U 1636-53. Their frequency 
separation decreases with the increase of one of the QPO frequencies. The black hole 
systems, on the other hand, e.g. GRO J1655-40, XTE J1550-564, GRS 1915-1-105, 
exhibit the kHz QPO pairs which seem to appear at a 3 : 2 ratio [li]. 

Interaction between the surface current density in the disk and the stellar magnetic 
field generating warps and subsequent precession instability in the inner accretion disk 
produces motions at low frequencies. This is similar to the Lense-Thirring precession 
|2]. This can explain the mHz QPOs for strongly magnetized neutron stars [i3J. The 
effect of nonlinear coupling between g-mode oscillations in a thin relativistic disk and 
warp has also been examined [0] for a static compact object. Recent observations ^ 
indicate a strong correlation between low and high frequency QPOs holding over mHz to 
kHz range which strongly supports the idea that QPOs are universal physical processes 
independent of the nature of the compact object. The correlation is also explained in 
terms of centrifugal barrier model [@]. Indeed earlier the correlation was shown in terms 
of the effective boundary wall created by the strong centrifugal force in the disk [3] . 

Accretion dynamics is a nonlinear hydrodynamic/magnetohydrodynamic phe- 
nomenon. It was already shown that QPOs may arise from nonlinear resonance 
phenomena in an accretion disk (e.g. [Hi @]) occurred due to resonance between 



epicyclic motions of accreting matter. However, the separation between vertical and 
radial epicyclic frequencies increases mostly as a function of either of the epicyclic 
frequencies contradicting the observed QPO feature. 

We aim at establishing a global model based on higher order nonlinear resonance 
theory to describe black hole and neutron star QPOs together and then try to predict the 
spin parameter/frequency of compact objects. The model predicts, alongwith important 
QPO features, the spin parameter of black holes as well as spin frequency of a neutron 
star. 

MODEL 

A system of degrees of freedom has N linear natural frequencies v,-; i = 1,2,...N fl^. 
These frequencies have commensurable relations which may cause the corresponding 
modes to be strongly coupled and yield an internal resonance. If the system is excited by 
an external frequency v*, then the commensurable relation exhibiting resonance might 
be 

N 

av* = biVi, with a + Y,\bi\ = j, (1) 

i=l i=l 

apart from all the primary and secondary resonance conditions cv* = dVm, where 
a, Z?;, c, d are integers and j = k+l when k is the order of nonlinearity. 

Now we consider an accretion disk with a possible higher order nonlinear resonance. 
The resonance is driven by the combination of the strong disturbance created by the 
rotation of the compact object with spin and the existent (weaker) disturbances in the 
disk at the frequency of the radial (v^) and vertical (Vj.) epicyclic oscillations given by 
(e.g. 

Vr = y ^/a^^4(V^^^, V, = y^r2-4av^ + 3a2, (2) 

where A = — 2Mr-\-a^, a is the specific angular momentum (spin parameter) of the 
compact object and is the orbital frequency of the disk particle given by ItzVq = ^ = 
l/(r3/2 + a). 

We describe the system schematically in Fig.[TJ It is composed of two oscillators due 
to radial and vertical epicyclic motion with different spring constants. The basic idea is 
that the mode corresponding to v^- will couple to the ones corresponding to and 
igniting new modes with frequencies v^^ i P^s LlQi], where p is a number e.g. L/2 or L, 
when L is an integer, if the effect is nonlinear or linear [[lol lsl] respectively. While L = 1 
corresponds to the dominant interaction, modes with other L are very weak to exhibit 
any observable effect. Now at a certain radius in the nonlinear regime where Vi/2 (or Vj 
in the linear regime) is close to — Vy ifiol [sl], and V5/2 (or v.,) is also coincidentally 
close to the frequency difference of any two newly excited modes, a resonance may 
occur which locks the frequency difference of two excited modes at V5/2 (or v^). 

As the neutron star has a magnetosphere coupled with the accretion disk, the mode 
with Vs can easily disturb the disk matter. A black hole, on the other hand, with magnetic 
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FIGURE 1. Cartoon diagram describing coupling of various modes in an accretion disk and corre- 
sponding nonlinear oscillator The oscillators describing by spring constant ki and k2 indicate respectively 
the coupling of spin frequency v* of the compact object with mass M to radial (x) epicyclic frequency and 
vertical (z) epicyclic frequency of a disk blob with mass m, where bi and b2 represent corresponding 
damping factors respectively. 



field connecting to the surrounding disk, may transfer energy and angular momentum 



to the disk by the variant Blandford-Znajek process [|12l] what was already verified by 
XMM-Newton observation (e.g. 1131]). This confirms possible ignition of new modes in 
a disk around a spinning a black hole. 

Therefore, rewriting ([T]) for an accretion disk we obtain 



n V5 m 

v. + -v, = v,-y + -v, (3) 



with —bi = b2 = 2 and m,n are integers. Now we propose the higher and lower QPO 
frequency of a pair respectively to be 

v/, = v,+ -v„ vi = Vz-—. (4) 

Hence, we understand from and ^ the order of nonlinearity in accretion disks 
exhibiting QPOs k = n — m + 4-. 

At an appropriate radius where — v,- ~ Vs/2 and resonance is supposed to take 
place, Av = — V/ ~ Vs/2 for n = m = I, and Av ~ v^. for n = m = 2. n = 1 implies a 
nonlinear coupling between the radial epicyclic mode and the disturbance due to spin of 
the neutron star, which results in Av locking in the nonlinear regime with m = 1. On the 
other hand, n = 2 implies a linear coupling resulting A v locking in the linear regime with 
m = 2. For a black hole, however, in absence of its magnetosphere, disturbance and then 
corresponding coupling may not be strongly nonlinear and occurs with the condition 
Vz — yr^ v.? resulting the resonance locking at the linear regime with n = m = 2, which 
produces Av ^ v., (sometimes ~ 2Vs/3). However, if we enforce the resonance strictly 
to occur at marginally stable circular orbit, then it occurs in the nonlinear regime with 



Av ^ V5/2 (sometime ~ Vs/5) and — V,- ^ v^- for n = m = l.ln principle, there may 
be possible resonances with other combination of n and m (e.g. n = 2, m = 1) which are 
expected to be weak to observe. 

Once we know the spin frequency v., of a neutron star from observed data we can 
determine specific angular momentum a (spin parameter) with the information of equa- 
torial radius R, spin frequency v,, mass M, radius of gyration Rq. If we consider the 
neutron star to be spherical in shape with equatorial radius R, then the moment of inertia 
and spin parameter are computed as / = MRq, a = IQ.sC / GM^ , where Q-s = iKVs, G is 
the Newton's gravitation constant and c is speed of light. We know that for a solid sphere 
Rq = 2R^/5 and for a hollow sphere Rq = 2R^/3. However, in practice for a neutron 
star Rq should be in between. Moreover, the shape of a very fast rotating neutron star 
is expected to be deviated from spherical to ellipsoidal. Hence, in our calculation we 
restrict 0.35 < [Rc/Rf < 0.5. 

On the other hand, for a black hole QPO a is the most natural quantity what we 
supply as an input. The corresponding angular frequency of a test particle at the radius 



of marginally stable circular orbit r,„s in spacetime around it is then given by 



g^t{r = r^s) 2a 
LIbh = 2;rv, = — — -— — - = 2 10 2 ' 

,[r — fms) ^ms'^nis^ -T Z,a 



where and light inside rms is practically not expected to reach us. 



RESULTS 
Neutron stars 

The choice of mass M — I.AMq and (Rc/R)^ = 0.4 does not suffice the observation 
mostly. Hence, we consider other values of parameters given in TABLE 1. 

For 4U 1636-53, our theory has an excellent agreement with observed data in accor- 



dance with realistic EOS [|l4ll shown in Fig. [2^ and given in TABLE 1. On the other 
hand, for 4U 1702-429, if the star is considered to be an ellipsoid and/or not to be a solid 
sphere and/or to have mass M < I .4-Mq, then our theory has an excellent agreement with 
observed data with realistic R [14] shown in Fig. [2b- Similarly for other neutron stars 
showing twin kHz QPOs, results from our model have good agreement with observed 
data (not discussed in the present paper but would be described elsewhere in future), 
which are the beyond scope to discuss in the present paper. 



Estimating spin ofSco X-1 

The spin frequency of Sco X-1 is not known yet. We compare in Fig. [2]: the observed 
variation of frequency separation as a function of lower QPO frequency with that 
obtained from our model. We find that mass of Sco X-1 must be less than 1.4M0 and 
results with sets of inputs with smaller and M fit the observed data better and argue 
that Sco X-1 is a slow rotator with ~ 280 - 300. 




FIGURE 2. Variation of QPO frequency difference in a pair as a function of lower QPO frequency for 
(a) 4U 1636-53, (b) 4U 1702-429, (c) Sco X-1. Results for parameter sets given in TABLE 1 from top to 
bottom row for a particular neutron star correspond to solid, dotted, dashed (4U 1702-429 and Sco X-1), 
long-dashed (Sco X-1) hues. The triangles are observed data points. 



TABLE 1. Vs is given in unit of Hz, M in M© , R in km, radial coordinate where 
QPO occurs, rQPo, in unit of Schwarzschild radius. 



neutron star 




M 




n.m 


R 


range of rgpo 


4U 1636-53 


581.75 


1.4 


0.5 


1 


16.5 


6.1-1.1 


4U 1636-53 


581.75 


1.1 


0.4 


1 


14.3 


8.2-9.3 


4U 1636-53 


581.75 


1.2 


0.35 


1 


16.8 


7.7-8.7 


4U 1702-429 


330.55 


1.0 


0.5 


2 


18.8 


10.2-11 


4U 1702-429 


330.55 


0.83 


0.35 


2 


18.5 


11.7-12.5 


estimated 


Sco X-1 


422.0 


0.9 


0.5 


1 


17.4 


9.4-12.5 


Sco X-1 


540.0 


1.2 


0.4 


1 


16.0 


8-9.1 


Sco X-1 


280.0 


0.8 


0.5 


2 


17.5 


11.1-14.4 


Sco X-1 


292.0 


0.81 


0.35 


2 


18.6 


11.2-14 



Black holes 

The possible mass or range of mass of several black holes is already predicted from 
observed data. However, the spin of them is still not well established. By supplying 
predicted mass and arbitrary values of a our theory reproduces observed QPOs for 
GRO J 1655-40 and GRS 1915+105 with their 3 : 2 ratio for n = m = 2 at a radius 
outside rms given in TABLE 2. However, if we enforce QPOs to produce at r^^ strictly, 
then they produce at a higher a fox n = m= 1 . Similarly, results from our model for other 
black holes showing twin kHz QPOs have good agreement with observed data (would 
be discussed elsewhere in future), which are the beyond scope to discuss in the present 
paper. 



TABLE 2. Vi^h are given in unit of Hz, M in Mq, tqpq and its distance from marginally stable orbit, 
Ar, are expressed in unit of Schwarzschild radius. 



black hole M 


a 






^QPO 


Ar n,m 


GRO J 1655-40 


estimated 


theory/observation 


theory/observation 






6-7 


0.737-0.778 


450/450 


300/300 


4.93-4.25 


1.71-1.23 2 


7.05 


0.95 


451.31/450 


299.04/300 


1.94 


1 


GRS 1915+105 


estimated 


theory/observation 


theory/observation 






10-20 


0.606-0.797 


168/168 


113/113 


7.38-3.9 


3.58-0.98 2 


18.4 


0.95 


167.35/168 


114.61/113 


1.94 


1 



SUMMARY 

We have prescribed a global QPO model based on nonlinear resonance mechanism 
in accretion disks. Based on this we have predicted the spin parameter/frequency of 
compact objects. The model has addressed, for the first time to best of our knowledge, 
the variation of QPO frequency separation in a pair as a function of the QPO frequency 
itself for neutron stars successfully. We argue that Sco X-1 is a slow rotator. 

We have addressed QPOs of black holes as well and predict their spin parameter {a) 
which is not well established yet. According to the present model, none of them is an 
extremally rotating black hole. As our model explains QPOs observed from several black 
holes and neutron stars including their specific properties, it favors the idea of QPOs to 
originate from a unique mechanism, independent of the nature of compact objects. 
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